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Abstract 

The set-up of the QCD Schrodinger functional (SF) on the lattice with staggered quarks requires 
an even number of points L/a in the spatial directions, while the Euclidean time extent of 
the lattice, T/a, must be odd. Identifying a unique renormalisation scale, L = T, is then 
only possible up to 0(a) lattice artefacts. In this article we study such lattices in the pure 
SU(3) gauge theory, where we can also compare to the standard set-up. We consider the SF 
coupling as obtained from the variation of an SU(3) Abelian and spatially constant background 
field. The 0(a) lattice artefacts can be cancelled by the existing 0(a) boundary counterterm. 
However, its coefficient, ct, differs at the tree- level from its standard value, so that one first needs 
to re-determine the induced background gauge field. The perturbative one-loop correction to 
the coupling allows to determine Ct to one-loop order. A few numerical simulations serve to 
demonstrate that residual cutoff effects in the step scaling function are small in both cases, 
T = L ziz a and comparable to the standard case with T = L. 



1 Introduction 



Renormalisation schemes based on the Schrodinger functional [l-3j have come to play an impor- 
tant role in lattice QCD and in Technicolor inspired models of electroweak symmetry breaking. 
The Schrodinger functional (SF) is used to define a technically convenient, intermediate renor- 
malisation scheme, where the scale is set by the finite space-time volume, and all fermion masses 
are taken to vanish. In SF schemes, the non-perturbative scale evolution of the running coupling 
or multiplicative renormalisation constants of composite operators can then be constructed re- 
cursively (see [IfSj for an introduction). Note that the scale evolution is eventually obtained in 
the continuum limit and thus universal. Hence, the choice of the lattice regularisation for this 
part of the calculation is merely a matter of practical considerations. In particular, in order to 
minimise computational costs it is advisable to use some variant of either Wilson or staggered 
quarks. 

The Schrodinger functional in lattice QCD was originally obtained for Wilson type quarks |:2j . 
Its formulation for staggered/Kogut-Susskind quarks has been initiated in |6] and further stud- 
ied in [7J. For theories with a multiple of four massless fermion flavourq^, staggered quarks 
constitute an attractive alternative to Wilson quarks: no tuning is required to recover the chiral 
limit, and numerical simulations for a given lattice size are computationally less demanding. In 
fact, the SF coupling has already been studied with iVf = 16 [8j and A^'f = 8, 12 [9j fermion 
flavours, and a four-flavour QCD study is in preparation |10H12j . A drawback, however, consists 
in the observation that staggered fermions require lattices with an even extent L/a in the spatial 
directions, whereas the time extent T/a must be odd |6|7j . As the renormalisation scale must 
be given by a unique scale, // = the ratio T/L must be fixed and is usually set to T/L = 1. 
With staggered quarks, this can only be achieved up to corrections of 0(a) and the question 
arises how the ensuing 0(a) effects can be eliminated. For the running coupling. Heller proposed 
to take an average of observables on lattices with T = L it a, and he showed that this procedure 
is consistent to one-loop order of perturbation theory. The same recipe was then also applied 
in numerical simulations ^8|9j. However, it is not obvious how this procedure can be generalised 
to fermionic correlation functions |13] . 

It is thus desirable to find an alternative to the averaging procedure. It is the purpose of the 
present paper to show how this is possible by slightly modifying the approach to the continuum 
limit. More precisely, the limit should be taken at fixed T' /L where T' is either set to T -|- a or 
to T — a. Requiring 0(a) improvement of the pure gauge theory leads to modifications of the 
boundary 0(a) counterterm proportional to Ct already at the tree-level. This in turn affects the 
equations of motion for the gauge field and thus changes the minimum action configuration used 
to define the coupling in the SF scheme. We determine the new background gauge field and 
then perform a one-loop computation to calculate the boundary 0(a) improvement of the SF 
to this order for the pure gauge theory. A numerical simulation has been carried out to check 
for the size of residual cutoff effects in the step-scaling function for the SF coupling. Although 
our motivation for this work originates in the fermionic sector of the SF, we here only discuss 
the necessary modifications to the SF in the pure gauge theory [1]. The details and calculations 
pertaining to the fermionic sector with staggered quarks will be presented elsewhere |10ll4j . 

^For staggered quarks we here use the traditional term "flavour" rather than "taste". 
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This paper is organised as follows: after a brief review of the continuum Schrodinger func- 
tional and the definition of the renormalised coupling, we discuss its lattice regularization and 

the origin of the constraint on lattice sizes with staggered fermions (Sect. 2). In Sect. 3, the 
background gauge field is found with Ct left as parameter, which is then determined such as to 
cancel 0(a) effects in the action. We then proceed with the one-loop calculation of the running 
coupling (Sect. 4), and determine ct to this order, as well as the size of remaining lattice artefacts 
at one-loop order. In Sect. 5 we discuss the results of our numerical simulation at an interme- 
diate value of the running coupling. Conclusions arc presented in Sect. 6. For future reference, 
an appendix contains two tables with the raw data of our perturbative one-loop calculation. 

2 A short reminder of the Schrodinger functional coupling 

In this section we review some known facts and definitions which will be relevant for this paper. 
In the continuum the Schrodinger functional is formally defined as the Euclidean path integral, 

Z\C,C'\ = J £>[A^,V^]e-^-"*[^'^'^l, (2.1) 

with the Euclidean action 



cont 



/ dxo / d^xCix), (2.2) 
Jo Jo 

£(x) = -^tr{F^,{x)Fi,,{x)} + ^{x){ji,D^ + m)^{x). (2.3) 

Here, go denotes the bare coupling constant, F^,y is the field tensor associated with the gauge 
field A^, 

Ff^u = d^A, - d,A^ + [A^, A,] , (2.4) 

and = + A^ + id^/ L denotes the covariant derivative acting on the quark fields (including 
a constant U(l) background field with = (1 — 5^q)0). The space-time manifold is taken to 
be a hypcr-cylinder. In the spatial directions periodic boundary conditions are imposed on all 
fields. At the time boundaries the conditions for the fermionic fields read 

P+V'Uo=o = = P_VUo=T, VSP_U=o = = ViP+Uo=r, (2.5) 

with the projectors P± = ^(1 it 70). For the gauge field one has 

^k\xo=o = Cj^, ^fc|i;o=r = C'fe, A; = 1,2,3, (2.6) 

with the boundary gauge field configurations and C'^. Of particular interest are spatially 
constant Abelian fields, which take the form, 

Cfe = ;^<^, C'k = ^4>', A: = 1,2, 3, (2.7) 
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V-o, (j)'i = -rj-TT, 



with traceless and diagonal N x A^- matrices (j) and cp' . For = 3 colours we follow [151 
parameterise the diagonal elements by, 

3' 

vi^-^y '^2=^('^+0 (2-8) 

where ry and are 2 real parameters. In the temporal gauge the field equations with these 
boundary conditions are solved by, 

Bo = 0, Bk = Ck + YiC'k-Ck), fc = 1,2,3. (2.9) 

which corresponds to a constant chromo-electric field, 

Go. = 5oi3, = = fc = l,2,3, (2.10) 

all chromo-magnetic components of the field tensor being zero. The background field B is, up to 
gauge transformations, uniquely determined by the gauge boundary fields since it corresponds 
to the absolute minimum of the gauge action, 

3 

where p = T/L is the aspect ratio. The effective action for the SF can thus be taken to be a 
function of the background field, 

T[B\ = -\iiZ[C',C], (2.12) 
which admits a weak coupling expansion, 

T[B] \to[B] + Tr[B] + 0{gl), (2.13) 

with the lowest order term given by the classical action, Tq[B] = g^ScontiB]. Setting T = L, one 
remains with a single external scale, L, and a renormalised coupling g{L) can now be defined 
through 



dr 

dr] 



k \ ^Kj-p- - iyv{L)\ , k ^ " 



r]=0 



g'^{L) j ' dr] 



= 12tt. (2.14) 

r/=0 



Here the ry-derivative serves to eliminate any background field independent terms in the effective 
action. Moreover, it implies that the coupling is defined by an SF correlation function rather 
than by the SF itself, which renders it measurable by numerical simulations. In practice the 
parameter v will be set to zero. However, its coefficient in Eq. (j2.14p . v, defines a further 
observable, which can also be measured at = 0, based on the relation. 



^(L) = --—{ — 
k dv\ dr] 



ri=0 



(2.15) 

u=0 
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In numerical simulations, the scale evolution for the coupling is constructed non-perturbatively 
by computing the continuum step-scaling function 

a(tx)=5'(2L)L=g2(i). (2.16) 

Prescribing a value for u = (fiL) implicitly fixes the scale L, and thus the only free parameter in 
pure Yang-Mills theory. The coupling at the scale 2L or v{L) are then fixed and can be obtained 
by taking the continuum limit of the corresponding observables on the lattice. We will later use 
the deviation from their respective continuum limits (t{u) and 

u:{u) = v{L)\^=-g2^L), (2.17) 

to obtain an impression of the cutoff effects for the different lattice regularisations considered 
here. 

The SF coupling is gauge invariant and can be non-perturbatively defined on the lattice. 
For small volumes it can be related perturbatively to the MS scheme. This relation is known to 
two-loop order both for the pure SU(3) gauge theory [16] and QCD [TTJ. Setting q = 1/L and 
a = /At:, the result for pure SU(3) gauge theory is 

"ms(9) = + cia^il) + C2a^iQ) + O(a^), (2.18) 

with [16], 

ci = 1.255621(2), C2 = c? -Fl. 197(10). (2.19) 
For the perturbative checks in this paper we will also need the one-loop result for the observable 

v{L) m, 

uj{u) = uJi+U2U + 0{u^), a;i = 0.0694603(1). (2.20) 
Finally, we recall that the step-scaling function to first non-trivial order is given by 

a{u) = u + 2boln{2) X u"^ + 0{u^), bo = U / (iir)'^ . (2.21) 

Here, bo is the universal one-loop coefficient of the SU(3) /3-function, which, for the SF coupling 
is known to 3- loop order |16j . 

3 Variants of the pure gauge Schrodinger functional on the lattice 

We here summarise the basic definitions and properties of the SF in the SU(A^) lattice gauge 
theory. The change at 0(a) motivated by staggered fermions has tree- level consequences which 
will be worked out in the remainder of this section. 

3.1 The standard lattice framework 

The lattice formulation of the SF for pure SU(A^) gauge theories can be obtained from the 
transfer matrix formalism yy. When written as a path integral it takes the form 

Z[C,C'] = [ DlU]e-^^^\ (3.1) 
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The pure gauge action is given by 



% p 

where the sum runs over all oriented plaquettes p of the lattice and U{p) denotes the parallel 
transporter around p. Assuming Abelian boundary gauge fields Ck and C^, the boundary 
conditions, Eq. (12. 6p translate to 

U{x,k)U,=o = W{jc,k), U{x,k)U,=T = W'{x,k), (3.3) 

with the boundary link variables, 

T^(x, k) = exp (aCfc), VF'(x, A;) = exp {aCQ. (3.4) 

Finally, the weight factor, w(p), is set to 

{^Cs{go) for spatial plaquettes p at xq = 0,T, 
ct{go) for time-like plaquettes p attached to the boundaries, (3.5) 
1 otherwise. 

In the continuum limit the plaquette terms multiplied by ct and Cg respectively reduce to the only 
gauge invariant local boundary operators of dimension 4, which are allowed by the symmetries 
of the SF (repeated spatial Lorentz indices are summed over), 

tr{FofcFofc}, tT{FkiFki}. (3.6) 

There are no other sources for lattice artefacts linear in a, and one may thus cancel 0(a) effects 
in any on-shell quantity by properly adjusting these coefficients. Moreover, while ct contributes 
to all observables, the counterterm proportional to Cg vanishes for spatially constant Abelian 
boundary gauge fields. In this case, 0(a) improvement can be achieved with the appropriate 
choice of Ct alone. In perturbation theory, one has 

Ct(5o) = cr + c;%o' + --- (3.7) 

In the standard set-up of the SF one sets L = T and tree- level 0(a) improvement is achieved by 
setting Ct = 1. With this choice the classical equations of motion which follow from the lattice 
action are equivalent to 

d*P{x,fi) - d*P{x,fi)^ - ^tr |d*P(x,^) - d*P{x,fi)^^ = 0, (3.8) 

where 

3 

d*P{x,f^) = ^{P^^(x) - Ul{x - au)P^^{x - av)U^{x - ai>)} , (3.9) 

and the plaquette field is defined by 

PA^) = Uf,{x)U^{x + afi)Ul{x + au)Ul{x). (3.10) 
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A solution V^{x) to Eq. (j3.8p which is unique up to gauge equivalence, is referred to as the lattice 
background gauge field. For spatially constant Abelian boundary fields, the lattice background 
field is obtained by simply exponentiating the continuum solution, Eq. (|2.9|) . 

V^,{x)=exp{aBf,{x)), // = 0,...,3. (3.11) 

Moreover, a mathematical proof establishes that this solution corresponds to an absolute mini- 
mum of the action, provided the condition 

TL/a^ > {N - l)7r2max{l, iV/16}, (3.12) 

is met |T]. For = 3 colours, this means that S[V] represents an absolute minimum for 
TL/o? > 27r^, implying L/a > 4 on lattices with L = T. The action at the minimum, 



(3.13) 



only differs at 0(a'^) from the continuum action, Eq. (j2.1ip . 

One should note, however, that these considerations depend on having set ct = 1. As we 
will now explain, staggered fermions require a modified approach to the continuum limit, which 
will lead to a different choice of ct. Consequently, the determination of the classical background 
field needs to be revisited. 



3.2 A constraint from staggered quarks 

A peculiarity in the formulation of the Schrodinger functional with staggered quarks has first 
been noticed in [6] and can also be traced back to the properties of the transfer matrix [18] : the 
lattice sizes needed in this case have even spatial extent L/a and odd temporal extent T/a. To 
understand this constraint, one needs to recall that the 4 flavours of Dirac (i.e. four-component) 
spinors are reconstructed from the staggered one-component fields living on the 2^ = 16 corners 
of the elementary hypercubes of the lattice. The reconstructed fields may be imagined to live on 
a coarse lattice with the doubled lattice spacing. The constraint arises from having to construct 
an integer multiple of four Dirac spinors from those staggered one-component fields which are 
integration variables in the functional integral (cf. the left panel of figure [1] for an illustration 
in 2 dimensions). Equivalently, the lattice geometry must be such that an integer number of 
hypercubes with dynamical field components is obtained. Note that the Dirichlet conditions for 
the one-component fields at xq = and xq = T naturally lead to a Dirichlet condition for half 
of the Dirac spinors living at the boundaries. This implicitly defines projectors in spin-flavour 
space, the form of which depends on the details of the reconstruction of the Dirac spinors. 
However, there is enough freedom to obtain four flavours of quarks satisfying the standard 
boundary conditions in terms of the projectors P± = ^(1 ± 70) [1917] . The reconstructed four- 
spinors thus live on a coarser lattice with size (L/2a)^ x (T/a — l)/2. 

A slightly less intuitive option to reconstruct the four-component spinors is illustrated in 
the right panel of figure [TJ It corresponds to an extension of the fine lattice by a time slice 
beyond each of the time boundaries. The part of the four-component spinors living on these 
added time-slices correspond to the non-Dirichlet components at the boundaries, which need not 
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be dynamical field component^. To treat this second case we also consider the coarser lattice 
with size {L/2a)^ x (T/a + l)/2, and, combining both options, we set 

T' = T + sa, s = ±l. (3.14) 

Note that in the absence of staggered quarks we may also set s = to recover the standard 
framework for the pure gauge theory. In the following we will thus try to take the continuum 
limit at fixed ratio, p = T' /L. In particular, for the computation of the SF coupling we set 
p = 1 exactly. This redefines the approach to the continuum limit, and we are thus led to 
discuss on-shell 0(a) improvement for this modified set-up. 





N on- dynamical 
fields 



Figure 1: The figure shows a 2-dimensional section of the SF with staggered fermions with 
L/a = A and T' /L = 1. The left panel shows the case T' = T — a (s = —1) the right panel 
corresponds to T' = T + a (s = 1). Thin lines represent the lattice on which the one-component 
staggered fermions live. The Dirac spinors are reconstructed from those one-component fields 
residing inside the circles and live on the effective lattice depicted by the thick lines. 



3.3 The equations of motion and the background field 

The background field is the field configuration which minimises the lattice action. Since the latter 
depends on ct, any change of this coefficient may modify the minimal action configuration. The 
equations of motion can easily be generalised to this case by simply including the weight factors 
Wfiuix) = w[P^u{x)] in the definition of the covariant divergence, Eq. ()3.9p . 



dl,P{x,p) = ^|u;^i,(x)P^i,(x) 

—w^u{x — ai')Ul{x — ai')P^u{x — az>)C/^(x — ai>)|. (3.15) 

The lattice action is stationary if and only if the traceless antihermitian part of d^P vanishes, 
i.e. the equations of motion read 

dlP{x, p.) - dlP{x, p)^ - Itr {dlPix, p) - dlPix, ^)t } = 0. (3.16) 



^This is in fact the situation in the SF with Wilson quarks, where the non-Dirichlet components at the 
boundaries are completely decoupled from the dynamical field components [2j. 
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In order to solve the equations we choose the temporal gauge, setting all temporal links to unity, 
Uo{x) = 1. Since the modification of the standard framework with T' = T amounts to an 0(a) 
effect, it seems plausible that the solution of the equations of motion will again be Abelian and 
spatially constant, at least for large enough lattice sizes. Hence, we make the ansatz 

V{x, k) = exp[a5fc(xo)], (3.17) 

and try to solve Eq. ()3.8p for a colour diagonal and spatially constant field i?fc(xo), subject to 
the boundary conditions 

5fc(0) = Ck = i(l)/L, Bk{T) = C'k = i<P'/L, (3.18) 

The coefficient Ct is left as a free parameter, which will be fixed later by demanding that the 
action is 0(a) improved. The field tensor associated with the plaquette field is defined by 

Pfj,^{x)\u-^v = exp[a^G^,,(xo)], (3.19) 
and is thus related to Bk{xQ) by the forward lattice derivative, 

G'ofc(xo) = doBkixo) = i [Bk{xo + a) - Bk{xo)] , (3.20) 

with all other components being zero. The boundary conditions are equal for all spatial indices 
k = 1,2, 3, so that we introduce the notation 

/(xo) = Gofc(xo). (3.21) 

Note that /{xq) is, for fixed xq, an anti-hermitian diagonal matrix in colour space with trace 
zero, i.e. its colour components /a(xo) are purely imaginary and must sum to zero. Eqs. (|3.16|) 
now reduce to T/a equations of the form, 

M(a;o) - ^tr [M(xo)] = 0, a<xo<T-a, (3.22) 

where the trace is over colour and the diagonal x A^-matrices M{xq) are given by 

{sinh [a^/(a)] — Ct sinh [a^/(0)] if xq = a, 

Ct sinh [a^f{T - a)] - sinh [a^fiT - 2a)] if xq = T - a, (3.23) 
sinh [a^/(2;o)] — sinh ^a'^f{xQ — a)] if a < xq < T — a. 

Expanding /(xq) in a power series in a, 

oo 

/M = E(f)"/(")(xo), (3.24) 



n=0 



and inserting this expansion in the equations of motion (j3.22p . one may show by induction that 
for all n, 

if a < xo < T - a 
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Hence one concludes that /(xq) is equal to the xo-independent constant matrix / except for a 
jump at the time boundaries by A/. To proceed further one needs to go back to the background 
field i?fc(xo) in the temporal gauge. One finds, 



Bk{0) 


= Ck 


Bk{xo) 




Bk{T) 





2 r 2 



if a < xq < T — a 



(3.26) 



and the jump in the field tensor at the boundaries is given by 

T 



A/ 



(3.27) 



2a 2a' 

Thus, for given boundary gauge fields and lattice geometry the background field and its field 
tensor are determined by the traceless and diagonal colour matrix /, which still needs to be 
computed. 

We now specialise to = 3 colours and insert the matrices / and A/ in Eq. (j3.22p . While 
M(xo) = for all o < xq < T — a, the 2 equations for M(a) and M(T — a) are identical and 
read, in terms of the colour components fa, 







1 



^ < ct sinh <^ o^//3 



1 + 



-Ct sinh < a^f, 



2 2a 



_ r_ 

~ 2^ 
+ ia- 



+ ia- 



2L 



smb. {a^f 13} 



2L 



+ sinh{a^/a} , 



(3.28) 



for a = 1, 2, 3. Taking into account that (/)[ - (p^ = -2((/)^ - cp^) = -2((/)3 - ^3) [cf. Eqs. (^Mj] 
one concludes that /2 and satisfy the same equation and hence /2 = fs- Tracelessness of / 
then implies 

/ = /2 xdiag(-2,l,l), (3.29) 

so that we are left with a single equation for /2. Introducing the real dimensionless variable (/? 
through 

/2 = i(p/L^, 

the equation to be solved for ip is 



Ct sin <{ 2j^K{ip) \ + Ct sin <{ 75 i^(9?) - sin <j 2^2 9^ 



a 



with 



K{ip) = if 



1 + 



pL 
2a 



+ 



L 
2a 



L2 



sin {j2^\=^, 



T'/L. 



(3.30) 
(3.31) 
(3.32) 



Next we expand ip in powers of a/L in analogy to Eq. (|3.24p and solve Eqs. (j3.3ip order by order 
in a/L. For the first 3 coefficients, we find. 



(0) 



(1) 



(2) 



{0)1, 



^2) 

kg 



-if ^ 



(2 + 5)ct - 2 

PCt 



(3.33) 

(3.34) 
(3.35) 
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Without further assumptions on ct or s there are corrections to all orders in a/L. We note, 
however, that for the standard case with s = and ct = 1, Eq. (13.3ip reduces to K{ip) = 99, 
which is exactly solved by 9? = ip^^\ In the general case one may obtain precise approximate 
solutions by pushing the Taylor expansion to higher orders. Alternatively, one may, for given 
numerical values of s, L/a, p, Ct and ij, find a numerical solution for (p (and thus for /) by 
applying the Newton procedure to Eq. (|3.3ip . 



3.4 Choice of ct 

The free parameter ct of the background field can now be determined by demanding that the 
classical action is 0(a) improved. This can be done analytically, as the dependence of the 
background field on Ct is known to 0(a) from Eq. (|3.34|) . First we insert the background field 
()3.26p into the action, and set T = T' — sa, 



S[V] 



-3[T' - (2 + s)a]L^ 



9'o 



6aL3 



Ct 



9i 



a=l 



y (/a + A/, 



Then, using the relations (|3.29|3.30p . and 



/ + A/ = -^K((^) xdiag(-2,l,l), 



the action can be written as a function of 99, viz. 

24L3 r 



S[V] 



2 2a 



sin ( j-^ip ) + 2 sin' 



(- 



Ct 



sin^ ( j^K{ip) ] +2 sin 



2L2 



Kiip) 



Then, with 



.(0) 



ip = if 

and Taylor expanding in a/L, one obtains 



(0) 



Ct 



l + ^ks + 0(a2) 



S[V] = ^ (^(°) 



l + -ks + 



(3.36) 



(3.37) 



(3.38) 



(3.39) 



0(a2)} = S,ont[B\ {1 + f + 0(a2)} . (3.40) 



In order to reduce the corrections to 0(0^) the coefficient kg, as defined in Eq. (j3.34p must 
vanish, which is achieved by setting 



Ct 



2 + s 



I if s = 1, 

1 if s = 0, 

2 if s 



(3.41) 



-1. 
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Remarkably, with this choice for ct, the action of the background field is not just 0(a) improved 
but cutoff effects only start at 0(a'^), in all three cases s = ±1, 0. In fact, the background field 
is only mildly distorted by lattice artefacts, as shown by the Taylor expansion, 

^k^.n.^s, = ^^'^ + (f )%^^^ + (f )%^^^ + (1)'%^"^ + O (a-) , (3.42) 
where (p^^^ is given in Eq. (I3.33P and 

^(^) = -M^.(. + 4)(. + 2), (3.43) 

^^'^ = - ^^l'^'^/ s{s + 4)(. + 2) (||.^ + f . + 17) , (3.44) 

= ^^^.^(. + 4)^(. + 2)^ = -|^(^(^))^ (3.45) 

The expansion to this order yields rather precise approximations to the full solution of Eq. (|3.3ip . 
For instance, with p = 1, we conclude that (p/(p^^^ is obtained with double precision provided 
L/a > 12 and L/a > 21 for s = — 1 and s = 1 respectively. 

In Section 4 we will also need the ?7-derivative of the background field, i.e. of (p. Given 
ip, this can be obtained by differentiating Eq. ()3.3ip with respect to rj. Alternatively one may 
obtain the a/L-expansion by differentiating the coefficients in Eqs. (j3.43ll3.45]) . Recalling that 
(p'2 — (p2 = ^ + T^/'^i the ?7-derivatives of the coefficients are given by 

= V> = ^. = J^. (3.46) 

3.5 A numerical check 

The above solution of the classical field equations was obtained using the hypothesis that the 
background field is Abelian and spatially constant. While this can be rigorously established in 
the standard set-up, we here rely on the plausibility of the assumption that an 0(a) change in 
the set-up can only have a small impact, at least for large enough lattice sizes. To check the 
hypothesis on small size lattices {L/a = 4 — 8) , we have tried to find the absolute minimum of the 
action numerically, by starting from a random gauge configuration. Sweeping through the lattice 
we minimised the plaquettes with respect to the SU(3) link variables. In order to avoid getting 
trapped in secondary minima of the action we also used over-relaxation steps which change the 
configuration whilst leaving the action invariant. The obtained results agreed numerically very 
well with the expectation from the Abelian background field. Furthermore, the background field 
itself was found to be of the expected form after transformation to the temporal gauge. We take 
this as evidence that the true minimal action configuration is indeed the Abelian and spatially 
constant solution described above. 

4 One-loop results 

The perturbation expansion is analogous to the calculation in SU(2) which has been described 
in the original paper [1] and we will be rather brief on the technical details. Given the one- 
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loop results we obtain the boundary improvement coefficients to this order and we can check 
remaining lattice artefacts in the step-scaling function for the coupling and in v{L). 



4.1 The running coupling at one- loop order 

So far our considerations have been purely classical. Dealing with quantum corrections the 
improvement coefficient ct will become a function of the bare coupling, with an expansion of the 
form, 

Ct(5o) = cl°)+5oM'^ + 0(ffo^). (4.1) 

Our classical considerations in the previous section correspond to setting go = i.e. Ct = c[^\ 
and c[^^ is thus given by Eq. (|3.41|) . The perturbative expansion of the effective action on the 
lattice takes the same form as in the continuum, Eq. (|2.13p . with the first 2 terms given by 



To[B] = <7o'^mL^=,(o) 



cf^ro;ct[^] -lndetAo + -lndetAi. 



(4.2) 
(4.3) 



Here, the fluctuation operators Aq and Ai appear in the Gaussian parts of the action in the 
Faddeev-Popov ghost and the gluon fields, respectively. The counterterm oc c[^^ is specified by 



oct <^t-'^t a-^ 



sm 



L2' 



K{ip) +2sin2 



2L2 



(4.4) 



where if is the solution of Eq. (I3.3ip with ct = cf^- Setting T' = L (i.e. p = 1) the definition of 
the SF coupling on the lattice 



dr,T[B] 



ri=v=0 



is such that g{L) = qq holds exactly at lowest order. Hence, the normalisation constant, 



a,ro[s] = 12 



2 



sm I — <^ 1+2 sin ( 2-pTip 



(4.5) 



(4.6) 



must be computed including the lattice artefacts and differs from its value in the continuum 
limit, k = 127r [cf. Eq. (l24l]l ]. by terms of O(a^). 

To one-loop order the perturbative relation to the bare coupling then reads. 



g\L)=g^^ + mi{L/a)-c\ 



(1) d.,T0;cAB] 



5o + 0(<7o'), 



r/=Oy 



d^To[B] 

where the coefficient mi{L/a) is a sum of two contributions, 

mi = [ho - ^/iij^^o- 
which derive from the Fadeev-Popov ghost and gluon fluctuation operators, 

dr, In det A 



d^ro[B] 



3 = 0, 1. 



(4.7) 



(4.9) 



r7=0 
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We have computed these contributions numericahy for lattice sizes ranging from L/a = 4 to 
L/a = 64, using two completely different techniques. One of us has performed the computation 
along the lines of refs. |l|2Uj . using recursion relations for the finite difference operators Aq 
and Ai, in order to compute their determinants in each sector of fixed spatial momentum and 
colour quantum numbers. In addition, an independent computation was carried out based on 
the automated perturbation theory described in refs. [21I22| . Perfect numerical agreement was 
found between the two methods, up to rounding errors. We also verified gauge invariance of 
mi{L/a) by computing it at different values of the gauge fixing parameter. For future reference 
the numerical results for mi{L/a) are collected in Tables [6] and [7] of Appendix A, for all three 
choices of s = ±1, 0.. To estimate the numerical precision we have compared the two calculations 
and find agreement for 15 digits at L/a = A which reduces to about 12 digits at L/a = 64, due 
to the accumulation of rounding errors. 

4.2 Determination of c[^^ 

From Symanzik's analysis |23|24j of the cutoff dependence of Feynman diagrams on the lattice, 
one expects mi{L/a) to have an asymptotic expansion of the form 

oo 

mi(L/a) Y.^a/Lr [r„ + s„ln(L/a)] . (4.10) 

n=0 

The series is logarithmically divergent, since mi (L/a) relates a renormalised to a bare coupling. 
For the coefficient of the divergence one expects sq = 26o, where 6o is the one- loop coefficient 
of the SU(3) /3-function, Eq. ()2.2ip . The divergence could be eliminated by renormalising the 
bare coupling e.g. by minimal subtraction of logarithms |25|26j . If one passes directly to the 
MS-scheme of dimensional regularisation, one finds that the coefficient ci in Eq. ()2.18p is related 
to ro in Eq. ()4.10p through 

ci = di(l)|jv=3 -47rro, (4.11) 
di{l) = -^ + kiN, A;i = 2.135730074078457(2), 

where the coefficient di{l) is taken from ref. [27|. Within the numerical precision we reproduce 
the known result for rg of ref. |15j in all three cases, s = ±1, 0. 

The coefficient ri multiplies cutoff effects of order a/L which ought to be cancelled by 
correctly choosing c^^-*. Hence we need to evaluate the coefficient of c|^'* in Eq. (|4.7|) to 0(a). 
To this order, the ct-dependence of Fq can be inferred from Eq. (j3.36p and is given by 

ro;cAB] = d,, [glS[V]]^^^^^, = glS,ont[B]j [dcM,.,'fi) + 0{a^), (4.12) 

and dcf^kg = 2/{pc^). This shows explicitly that this term is of 0(a). Performing the ry-derivative 
and requiring the absence of 0(a) terms in the one- loop relation between the couplings, leads 
to the condition. 

Note that the coefficient si of the term {a/L)x In (L/a) is expected to vanish, provided tree-level 
improvement is correctly implemented. Following the method presented in the appendix of [17j, 
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the first few coefficients were extracted numerically. The outcome of this analysis for the cases 
s = 0, ±1 is shown in Table [H In this table, coefficients shown with no errors have been assumed, 
entries x have also been fitted but are not listed here and the terms ~ are included one at a time 
for the error analysis. As can be seen from the results, all expectations are confirmed within 





So 


ro 




n 


S2 


r2 


S3 




s = -l 


0.13931(4) 


0.3683(3) 


-0.001(15) 


-0.23(7) 


X 


X 






22/(47r)^ 


0.368283(3) 


-0.0004(10) 


-0.230(6) 


X 


X 






22/(47r)^ 


0.3682818(7) 





-0.2318(3) 


X 


X 






s = 


0.13931(3) 


0.3683(2) 


-0.001(13) 


-0.17(7) 


X 


X 






22/(47r)^ 


0.368283(2) 


-0.0003(9) 


-0.176(6) 


X 


X 






22/(47r)^ 


0.3682817(7) 





-0.1779(3) 


X 


X 






s = 1 


0.13931(2) 


0.3683(1) 


0.001(7) 


0.12(3) 


X 


X 






22/(47r)^ 


0.368280(1) 


0.0007(5) 


0.120(3) 


X 


X 






22/(47r)^ 


0.368283(1) 





0.1232(4) 


X 


X 







Table 1: Asymptotic expansion coefficients in mi{L/a), for s = ±1,0. 



the errors. In particular we confirm the known result for the case s = [15], and the universal 
coefficients rg and sq agree for all choices of the parameter s. As expected, si is found to be 
compatible with zero. Applying formula (I4.13P with p = 1 we obtain 

0.0274(1) if s = 1, 

-0.08895(15) if s = 0, (4.14) 
-0.4636(6) if s = -1, 

The result for s = agrees within errors with the more precise value cj^^ = -0.08900(5) quoted 
in ref. [15]. Note that ct is the coefficient of a local counterterm and thus cannot depend on 
global space-time properties such as the aspect ratio p. As a further check we also produced some 
data with p = 1/2. Within the numerical accuracy, the coefficient ri is found to be proportional 
to 1/p, which cancels the explicit factor p in Eq. (j4.13p . and renders c[^^ p-independent, as 
expected. 



* 2 \2 + sJ 



4.3 Residual cutoff effects in the step scaling function 

In the previous subsection we have removed the 0(a) boundary lattice artefacts at one-loop 
order in perturbation theory. Nevertheless, higher order lattice artefacts are still present. We 
quantify these by studying the relative deviation of the lattice step scaling function, Y,{u,a/L) 
from its continuum counterpart, cr{u) [cf. Eq. (j2.2ip ]. 

S{u,a/L) = ^(^' Q/-^) - ^(^) ^ Si{a/L)u + 62{a/L)u'^ + O(-u^). (4.15) 
a{u) 



Defining 

mi(L/a) = mi(L/a) - c^^ ^^^^^ 



(4.16) 

ri=0 
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a one-loop computation yields 



8i(a/L) =mi{2L/a) -mi{L/a) - 26oln(2), 



(4.17) 



which we here use to monitor higher order lattice artefacts. In practice, 5i{a/L) can also be 
used to cancel all lattice artefacts to O(n^) in non-perturbative estimates of the step-scaling 
function. However, this implies that the ct-counterterm must be evaluated exactly and not just 
in its asymptotic form as in Eq. (j4.12p . Using the field equation (j3.3ip . the resulting formulae 
are surprisingly simple, and one obtains, 



SO) L 



1 H P — On<4^ H 

2 ^2a 2a 



a{2 + s) 
L 2^" 



(4.18) 



In the case s = one has drj^p = 1/p from Eq. (j3.33p and the exact ct-counterterm is thus found 
to coincide with its asymptotic a/L-term. For s = ±1 the numerical difference is very small, due 
to the mild cutoff effects in (/?, Eq. p.42p . and hence in dr^(p. In Table [2] and in Fig. [21 we show 
two sets of data for 6i{a/L), obtained by either setting c[^^ = or by choosing the correct value 
for c[^^ from Eq. ()4.14p . Asymptotically one expects 6i{a/L) to approach zero with a rate (x a/L 
and oc (a/L)'^, respectively. This is indeed observed in all cases. We also note that the data 
for all choices of the parameter s seem to be behave similarly, i.e. there is no clearly superior 
choice. 





<^i(a/^)l4i)=o 


Siia/L) 


L/a 


s = 1 


s = 


s = -1 


s = 1 


s = 


s = -1 


4 


-0.02036 


0.01192 


0.01575 


-0.00489 


-0.01033 


-0.01302 


6 


-0.01224 


0.01089 


0.01519 


-0.00195 


-0.00394 


-0.00409 


8 


-0.00876 


0.00918 


0.01252 


-0.00105 


-0.00194 


-0.00196 


10 


-0.00686 


0.00773 


0.01041 


-0.00069 


-0.00117 


-0.00117 


12 


-0.00564 


0.00663 


0.00887 


-0.00050 


-0.00079 


-0.00079 


14 


-0.00479 


0.00579 


0.00771 


-0.00039 


-0.00057 


-0.00057 


16 


-0.00416 


0.00513 


0.00682 


-0.00031 


-0.00043 


-0.00043 


18 


-0.00368 


0.00461 


0.00610 


-0.00025 


-0.00034 


-0.00033 


20 


-0.00329 


0.00418 


0.00553 


-0.00021 


-0.00027 


-0.00027 


22 


-0.00298 


0.00382 


0.00505 


-0.00018 


-0.00022 


-0.00022 


24 


-0.00272 


0.00352 


0.00464 


-0.00015 


-0.00019 


-0.00018 


26 


-0.00250 


0.00326 


0.00430 


-0.00013 


-0.00016 


-0.00016 


28 


-0.00232 


0.00304 


0.00400 


-0.00012 


-0.00013 


-0.00013 


30 


-0.00216 


0.00285 


0.00375 


-0.00010 


-0.00012 


-0.00012 


32 


-0.00202 


0.00268 


0.00352 


-0.00009 


-0.00010 


-0.00010 



Table 2: One-loop coefficients 5i{a/L) for s = 1,0, —1, without and with the one-loop countert- 

-0.089 [15]. 



erm oc c|^\ For s = ±1 we have used Eq. ()4.14p . for s = we set c^^ 
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s=-1 



0.02r 

0.015 
0.01 



-2- 

«r - 

-0.005- 
-0.01 - 
-0.015^ 

0.015^ 
0.01 - 
0.005- 

0- 

-0.005- 
-0.01 - 
-0.015^ 

0.005^ 

0- 

-0.005- 
«, -0.01 - 
-0.015- 
-0.02- 



* Unimproved o Improved 



0.005- ^ 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 

{a/lf 

s=0 



* Unimproved o Improved 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 

(a/Lf 

s=1 



* Unimproved o improved 



30 o o 



0.01 0.02 0.03 0.04 0.05 0.06 0.07 

(a/Lf 

Figure 2: Discretisation effects in the step-scaling function for s = —1,0, 1 from top to bottom, 
respectively. Stars show the result obtained before cancellation of the 0(a) part by the pure 
gauge boundary counterterm, c[^'^ and circles the results after the cancellation. 
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4.4 Residual cutoff effects in v{L) 

At fixed u = g'^{L) one may also study the lattice artefacts in the observable v{L) introduced 
in Section 2. On the lattice we define, 

n{u,a/L) = v{L)\^=-g2(^L) =ni{a/L) + n2{a/L)u + 0{u^), (4.19) 

where Q{u, a/L) converges to the universal function u}{u), Eq. (j2.17p in the continuum limit. In 
order to compute the one- loop coefficient ^}i{a/L) it is convenient to consider a finite difference in 
z^- values, rather than using Eq. (j2.15p . Both are equivalent since the i^-dependence is completely 
explicit in Eq. ()2.14p . Hence, 



v{L) 



1 / dT[B] 



_ dT[B] 



(4.20) 



ri=0,u=U2 



{1/2 - z^i)5^ro[-B]|^=o \ df] 

Choosing i^i = and = 1 one obtains 

ni{a/L) = m'{=^{L/a) - mi{L/a), (4.21) 

with the one- loop coefficient [cf. Eq. (j4.9p ]. 

mr'=[ho-hhi],=,- (4.22) 

Note that the Ct-counterterm is i^- independent and does therefore not contribute in Eq. (|4.2ip . 
We have computed m^'^^ for even lattice sizes up to L/a = 64. By combining with the u = 
data we have checked that Qi{a/L) converges to the universal continuum value, ui of Eq. ()2.20p 
for all cases s = ±1,0. 

To study the cutoff effects we follow ref. \T5] and define, 

a/L) ^ ^M^-^iu) ^ ^ e2{a/L)u + Oiu"). (4.23) 

0J[U) 

At one-loop order we obtain 

ei(a/L) = 0(a/L)/a;i - 1, (4.24) 

and we collect our numerical results in table [3l The continuum limit is always approached with 
a rate of O(a^), as expected. While the cutoff effects in ili are sizeable, we do not observe a 
striking difference between the cases s = 0, ±1. 



5 A non-perturbative check 

In order to check our set-up non perturbatively, we have carried out numerical simulation to 
compute the step scaling function and the observable v. To generate a representative ensemble 
of gauge configurations, we used an algorithm consisting in a combination of local heat bath 
sweeps (HB) with a number of over-relaxation sweeps (OR). Measurements are performed for 
every cycle, where a cycle is defined by, 

1 cycle = 1 HB + iVoR x OR, (5.1) 
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L/a 


s = -1 


s = 


s = 1 


4 


0.21001 


0.18806 


0.10209 


6 


0.06874 


0.06773 


f~\ r\ AAA i~\ 

0.04449 


8 


0.03365 


0.03360 


0.02419 


10 


0.02040 


0.02039 


0.01562 


12 


0.01380 


0.01380 


0.01104 


14 


0.00999 


0.00999 


0.00825 


16 


0.00758 


0.00758 


0.00641 


18 


0.00595 


0.00595 


0.00513 


20 


0.00480 


r\ r\r\ A or\ 

0.00480 


0.00420 


22 


0.00395 


0.00395 


0.00350 


24 


0.00331 


0.00331 


0.00297 


26 


0.00282 


0.00282 


0.00254 


28 


0.00243 


0.00243 


0.00221 


30 


0.00211 


0.00212 


0.00193 


32 


0.00185 


0.00185 


0.00171 



Table 3: One-loop coefficients ei[a/L) for s = —1,0, 1 

with A'^oR = L/{2a). To study the evolution of the running coupling, pairs of lattices with sizes 
L/a and 2L/a have been simulated. The results are shown in Table HI In order to extract the 
mean value of the observables, as well as their uncertainties and their integrated autocorrelation 
times, we have used the procedure advocated in ref. [28] . 

In Table [5] we show the results for the step scaling function at fixed u = 2.1. Figure [3] shows 
the cutoff dependence of the step scaling function at u = 2.1 and the continuum extrapolation. 
The results for s = from ref. |15j are shown for comparison. 



s 




L/a 


f 


Tint 


V 


Tint 


stat. 


+1 


7.2095 


6 


2.0966(31) 


0.64(1) 


0.0579(15) 


1.49(4) 


4700 X 16 


12 


2.7191(96) 


0.85(2) 


0.0501(24) 


1.89(8) 


4500 X 16 


-1 


7.1214 


6 


2.0929(45) 


0.58(1) 


0.0610(16) 


0.91(2) 


29700 X 4 


12 


2.739(15) 


0.64(2) 


0.0470(29) 


0.98(4) 


4700 X 16 


+1 


7.4218 


8 


2.0989(37) 


0.68(1) 


0.0602(17) 


1.69(6) 


5500 X 16 


16 


2.725(17) 


0.93(4) 


0.0524(37) 


1.74(10) 


9500 X 4 


-1 


7.3632 


8 


2.0968(47) 


0.58(1) 


0.0587(16) 


0.93(2) 


39700 X 4 


16 


2.747(18) 


0.59(2) 


0.0557(36) 


0.95(4) 


16700 X 4 


+1 


7.7447 


12 


2.0937(44) 


0.77(2) 


0.0578(20) 


1.90(6) 


7500 X 6 


24 


2.715(15) 


0.95(3) 


0.0502(34) 


1.90(9) 


19500 X 4 


-1 


7.6985 


12 


2.1017(55) 


0.61(1) 


0.0579(19) 


0.99(2) 


49700 X 4 


24 


2.772(19) 


0.67(2) 


0.0452(34) 


0.97(4) 


29700 X 4 



Table 4: Simulation parameters and results for and v, and their autocorrelation time Tint. 
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C\J 

II 




2.68 



0.05 



0.1 

a/L 



0.15 



0.2 



Figure 3: Cutoff dependence of the step scaling function at n = 2.1 and continuum extrapolation. 
For comparison, results for s = taken from ref. [15J are also shown. 



6 Conclusions 

Motivated by applications to lattice theories with staggered fermions, we have defined an 0(a) 
modified lattice set-up for the Schrodinger functional in the SU(3) gauge theory. Both perturba- 
tively and non-perturbatively the cutoff effects in the step-scaling function for the coupling are 
comparable to the standard set-up. Note that our approach can be applied to the SU(2) gauge 
theory with minor changes. However, for SU(A^) with > 4 this may or may not be the case, 
depending on the choices made for the boundary gauge fields. More precisely, assuming that 
the induced lattice background field remains Abelian and spatially constant, one may need to 
solve coupled equations in two or more variables replacing the single equation for y?, Eq. ()3.3ip 
of the SU(3) theory. 

To extend this framework to four-flavour QCD and other QCD-like theories with staggered 
fermions, one needs to discuss 0(a) boundary improvement with staggered fermions |13|14j . 
While our motivation originates in applications using staggered fermions, we note that there 
are other potential applications: for instance, the SF with Wilson or Ginsparg- Wilson fermions 
and chirally rotated boundary conditions is, for technical convenience, constructed with an 0(a) 
offset in the orbifold construction |30|31j . This entails tree- level 0(a) artefacts in the fermionic 
propagator which could be avoided by keeping T'/L fixed when taking the continuum limit, 
with T' = T ziz a. In any case, the new framework increases the flexibility for applications of 
the SF scheme. Furthermore, it may help to improve control over the continuum limit, e.g. by 
performing constrained continuum extrapolations of data for step-scaling functions, computed at 
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L/a 


s 


<j{u = 2.1) 


6 


+1 


2.725(11) 


D 


-1 

— 1 


2.751(17) 


8 


+1 


2.727(18) 


8 


-1 


2.752(21) 


12 


+1 


2.726(17) 


12 


-1 


2.769(21) 


cont. 


+1 


2.728(35) 


cont. 


-1 


2.783(44) 



Table 5: Continuum extrapolated results for the step scaling function at u = 2.1 . Central values 
and error are corrected by following the strategy in section 3.1 of ref. |29] . 

different values of the parameter s. This has already been attempted with staggered fermions 
[T2] and it can be applied more generally. 
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A Raw one-loop data 



LI a 

1 


s 


= 1 


s 


= 


s 


= -1 


4 





59867564635006385 





52976131887602069 





51922066331094546 


6 





64113765437070799 





59321232599981455 





58442701843014381 


8 





67488327354469054 





63824441829314597 





63153838271523941 


10 





70238903864308802 





67280038634701944 





66742121983892173 


12 





72546473348301643 





70067317431899355 





69618664881505825 


14 





74530462655573632 





72398165609402260 





72013453722120221 


16 





76269338777455524 





74399271936483262 





74062578372376441 


18 





77816541775427566 





76151561542170611 





75852242009405113 


20 





79209901867921951 





77709645174594405 





77440237357962767 


22 





80477118528362614 





79112035275707051 





78867107103739184 


24 





81639058776392414 





80386875084304289 





80162350146808351 


26 





82711823835863092 





81555339356121976 





81348080743360567 


28 





83708093056087357 





82633751751530950 





82441294090339149 


30 





84638026935694360 





83634957239642928 





83455327855159342 


32 





85509891481710033 





84569242491756803 





84400838361422786 


34 





86330500808698590 





85444971446735487 





85286472295855835 


36 





87105537841668160 





86269035439316518 





86119340950763814 


38 





87839791234295844 





87047179120264344 





86905362653451155 


40 





88537333426399287 





87784241080234770 





87649514951940602 


42 





89201656536325291 





88484334593899240 





88356023620527862 


44 





89835777515201167 





89150985494698484 





89028506544613721 


46 





90442320537524014 





89787238814607281 





89670084803647420 


48 





91023582290682148 





90395742302994047 





90283469523001835 


50 





91581584250088787 





90978812583979161 





90871030564770146 


52 





92118114931205926 





91538488106299742 





91434851426796011 


54 





92634764338015540 





92076571925410400 





91976773540748381 


56 





93132952274785759 





92594666572356148 





92498432332375374 


58 





93613951787157274 





93094202701724371 





93001286814371605 


60 





94078908704027068 





93576462803843164 





93486644054503752 


62 





94528858033171669 





94042600966826691 





93955679547364682 


64 





94964737799172085 





94493659452353015 





94409454285563946 



Table 6: The one-loop coefficient, mi{L/a), Eq. (14. Sp for s = 1,0, —1. By comparing results of 
two independent calculations we estimate that the number of significant digits decreases due to 
rounding errors from 15 for L/a = 4 to about 12 for L/a = 64. 
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L/a 


s 


= 1 


s 


= 


s 


= -1 


5 





62117136258700821 





56468083125839683 





55461713132090261 


7 





65893111142441375 





61737811787009622 





60974739393470162 


9 





68928647150162301 





65654153428089878 





65057003087364097 


11 





71439398959413626 





68741301652316919 





68252028857956904 


13 





73573276395961912 





71280479111934452 





70866242179288317 


15 





75426748654695955 





73434096296455847 





73074988732724165 


17 





77064247215528033 





75302638176979129 





74985728431314708 


19 





78530529883645214 





76952176876897567 





76668599155622038 


21 





79857842470262715 





78428347860819714 





78171762080518526 


23 





81070148601313597 





79763941632436668 





79529658268786343 


25 





82185727793917046 





80983289063631291 





80767742403345886 


27 





83218834749448024 





82104930186717609 





81905346010774183 


29 





84180796784671203 





83143311025949667 





82957488596710694 


31 





85080762207365805 





84109903085829275 





83936067325760364 


33 





85926224466335029 





85013965596714345 





84850663994188621 


35 





86723397937992463 





85863078810432712 





85709107745685024 


37 





87477492945292666 





86663526032990589 





86517876996699599 


39 





88192920731676446 





87420573028207099 





87282392619177697 


41 





88873448731299464 





88138676139487858 





88007235812526203 


43 





89522319912792745 





88821639863152864 





88696312708410143 


45 





90142345720374637 





89472737905126752 





89352980585123832 


47 





90735979318226442 





90094807414573837 





89980145942895982 


49 





91305373939615786 





90690323216425240 





90580341637673730 


51 





91852429830862484 





91261456924592534 





91155788213294810 


53 





92378832362922261 





91810124483365024 





91708443159609367 


55 





92886083231152201 





92338024750431470 





92240040837837180 


57 





93375526193982340 





92846671072289073 





92752125115677514 


59 





93848368458072277 





93337417324448625 





93246076251996867 


61 





94305698564346882 





93811479540720059 





93723133204075556 


63 





94748501439735390 





94269953998217965 





94184412261958892 


65 





95177671136930530 





94713832432409256 





94630922711111939 



Table 7: Same as table 6 for lattices with odd L/a. 
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